Fluctuations of a homeotropically aligned nematic liquid crystal in the presence of an 

applied voltage 
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We determined the refractive-index structure-factor S„(k) from shadowgraphs of fluctuations in 
a layer of a homeotropically aligned nematic liquid crystal with negative dielectric anisotropy in the 
presence of an ac voltage of amplitude Vo applied orthogonal to the layer. S'„(k) had rotational 
symmetry. Its integral P(Vo) and amplitude B(Vo) increased smoothly through the Freedericksz 
transition at Vo — Vf- Its inverse width £(Vo) and its relaxation rate Vo(Vo) had cusps but remained 
finite at Vf- The results are inconsistent with the critical mode at a second-order phase transition. 



A nematic liquid crystal (NLC) consists of elongated 
molecules that, because of their shape, tend to align lo- 
cally relative to each other. [l| The alignment direction 
is called the director n. When a NLC is confined be- 
tween parallel glass plated with a small spacing d between 
them, h is influenced by the interaction of the molecules 
with the glass surfaces. It is possible to prepare the sur- 
faces in such a fashion as to cause more or less uniform 
director alignment throughout a thin sample. When n 
is orthogonal (parallel) to the surfaces, the alignment is 
called homeotropic (planar). We note that planar align- 
ment determines a preferred direction in the plane of the 
sample and thus breaks the rotational invariance charac- 
teristic of an isotropic fluid. However, in the homeotropic 
case, which is of interest here, this rotational symmetry 
is preserved. 

Most physical properties of a NLC are anisotropic. We 
applied a voltage V = V2Vq cos(wi) between transparent 
indium-tin oxide (ITO) electrodes on the inner surfaces 
of the confining plates, and thus are interested in the 
anisotropy of the dielectric constant e a = eii — e_i_ where 
eu and e± are the dielectric constants parallel and perpen- 
dicular to the director. When e a < 0, the homeotropic 
state becomes unstable when Vo exceeds a threshold value 
Vf- The transition is known as the Freedericksz transi- 
tion (FT). @ Above Vf the director, in the ideal case, 
remains homeotropically anchored at the surfaces, but 
away from the confining plates it acquires a component 
n x in the plane (we arbitrarily choose the x-axis in the 
direction of this deformation). 

A stability analysis at the mean-field level of the FT 
was given by Hertrich et al. [3| They found that the tran- 
sition is of second order and that it occurs first for wave 
number k = 0. Here we are concerned with the time- 
dependent fluctuations of the director field [100 that 
are induced by thermal noise and not seen in global or 
time-averaged measurements 0. A spatial variation of 
the component of h in the plane of the sample is associ- 
ated with a variation of the vertical average of the refrac- 
tive index n and thus can be seen by the shadowgraph 
method. Q For the purpose at hand the method has 
the advantage that it sees only the fluctuations and not 



the much larger spatially uniform change of the refrac- 
tive index above the transition that influences the total 
transmitted intensity when the sample is viewed between 
crossed polarizers. |7j This is so because its sensitivity at 
small k is proportional to k 2 and thus vanishes for k = 0. 
From shadowgraph images we derived the structure fac- 
tor 5„(k, Vo) of the refractive- index field using the optical 
transfer function derived from physical optics. [lOj 

We looked for, but did not find, any hysteresis at 
Vf- For fluctuations near the critical point of an equi- 
librium system |ll| or a supercritical bifurcation of a 
non-equilibrium system [l2l Il3l IhH HH H(| the fluctu- 
ation power P of the critical mode should pass through 
a sharp maximum as a control parameter e, in our case 
equal to Vq /Vp — 1, increases from negative to positive 
values. Such a maximum was observed for instance near 
the onset of RBC.^3| We found that P varied smoothly 
and monotonically through e = 0, continuing to increase 
with e above the transition. The correlation length £(Vo), 
equal to the inverse width at half height of the azimuthal 
average S n (k) of 5 n (k), is expected to diverge at e = 0. 
We found that it had a sharp but finite cusp, indicat- 
ing a well defined transition with "rounding" confined 
to |e| < 0.01. The relaxation rate To(Vo) is expected to 
vanish at Vf] it was found to have a sharp down-pointing 
cusp at Vf but remained non-zero. The height B of the 
fluctuation spectrum should be pro port ional to the sus- 
ceptibility of the fluctuating mode and is expected 
to diverge. Nearly in proportion to P, it also varied 
smoothly and monotonically through Vp. We are forced 
to one or the other of two not very palatable conclusions. 
One of them could be that the observed fluctuations, al- 
though influenced by the transition, have an origin that 
is different from that of the critical Freedericksz mode. 
We find this unlikely because the shadowgraph method 
should reveal an average, orthogonal to the plane, of in- 
plane director fluctuations that are believed to be asso- 
ciated with the critical mode, and because the time scale 
of the fluctuations is comparable to the relevant director 
relaxation-time ~ 1 s. Alternatively, the measure- 
ments are consistent with a first-order phase transition 
in the presence of noise sufficiently strong to eliminate 
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hysteresis and to cause the transition to occur at the 
thermodynamic transition point where the free energies 
of the two phases are equal rather than at the point of 
absolute instability of the ground state. 

Our results differ form light scattering measurements 
on related but different systems that also have a FT. In 
a sample with planar alignment and e a > 0, Q and in a 
homeotropic sample in the presence of a magnetic field 
in the plane, [§) r (Vo) was found to vanish at Vf- These 
systems differed from ours in that the ground state below 
the FT was not rotationally invariant, either because of 
the director field or because of the magnetic field []| . 

The rotational invariance of our system below Vf is a 
feature in common with micro-crystallization of di-block 
co-polymers J 171 an d the onset of Rayleigh-Benard con- 
vection (RBC) ]H E3- These two belong to the Bra- 
zovskii universality class. Members of this class are 
expected to exhibit a fluctuation-induced first-order tran- 
sition even though at the mean-field level the transition 
is of second order. What distinguishes the FT from RBC 
and the co-polymer case is that it is an instability at 
k = 0. In the other systems the rotational invariance is 
broken at the transition by the selection of a direction in 
the plane due to the formation of a striped phase with 
a wave director of finite length; in the present case a 
spatially uniform domain, corresponding to k = 0, forms 
above Vf and the symmetry breaking is attributed to the 
existence of the nematic director within this domain. 

Regardless of the nature of the transition, for macro- 
scopic samples fluctuation effects should affect the nature 
of the phase transition only very close to the transition. 
One can estimate rou ghly that the relevant thermal noise 
intensity is given by [lj[ F t h — ksT/k^d. Here fc 33 ~ 
8.6 x 10~ 12 N is the bend elastic constant. Q For our 
d = 2.7 x 10 -5 m one has F th ~ 2 x 10~ 5 . The expected 
critical region has a width of order e c ~ F 2 / 3, ~ 7 x 10~ 4 . 
|l2l | If we interpret our results as a fluctuation-induced 
first-order transition, the data suggest e c = O(0. 1) and 
would require a noise intensity about three orders of mag- 
nitude larger than the theoretical estimate. We find it 
unlikely that there is an experimental noise source that 
can account for the observations because it would have 
to be additive (rather than multiplicative) and it would 
have to have the appropriate temporal and spatial spec- 
trum. 

The main sample was N-(4-methoxybenzylidine)-4- 
butylaniline (MBBA) with d — 27 /im, doped with 
less than 0.01% of tetra-n-butyl-ammonium bromide 
(TBAB), at 25°C driven at / = 100 Hz. The same re- 
sults within our resolution were obtained at / = 1000 
Hz. The sample was confined between glass plates 
covered on the inside by ITO electrodes that in turn 
were coated with a thin film of dimethyloctadecyl[3- 
(trimethoxysilyl)-propyl] ammonium chloride (DMO AP) 
to produce homeotropic alignment. The sample had a 
conductance a = 5 x 10 -7 (ohm m) _1 , but similar results 




FIG. 1: Left column: shadowgraph snapshots of an area of 
1.03 x 1.03 mm 2 . Right column: corresponding structure 
factors, averaged over 256 images, (a): Vo = 3.170 , (b): 
V = 3.256 Volt. 

were obtained with undoped MBBA and a = 3 x 10 -8 
(ohm m) . The apparatus, although different in detail, 
was equivalent to one described before. 0, Mea- 
surements between crossed polarizers for Vq > Vf showed 
that the imaged area of 1.03 x 1.03 mm 2 was within a sin- 
gle domain. Time sequences of Images ii(x, Vq,t), i — 
1, N, N — 256 (x = (x, y) is the two-dimensional po- 
sition vector in the plane of the sample) were acquired at 
intervals of three to fifteen seconds with various camera 
exposure times 0.025 < t < 8 s at numerous values of Vb 
over the range 2.5 < Vq ?S 3.6 Volt. No polarizers were 
used. At each Vo and r each Zj was divided, pixel by 
pixel, by the average /o( x , Vo, r) of Ii. This removed any 
time-independent structure but left most of the fluctua- 
tions. The shadowgraph signal X = Ii/ Iq — 1 was Fourier 
transformed, and the structure factor S(k, Vo,t), equal 
to the square of the modulus of the Fourier transform, 
was obtained. 

Figure^shows grey-scale renderings, rescaled by their 
own variances, of X (left column) and S (right column). 
Here (a) is well below and (b) is about one percent above 
V F = 3.22 Volt (e ~ 0.02). There is no obvious change 
of X and S with Vq. It is remarkable that S is rota- 
tionally symmetric even above Vf where the formation 
of a Freedericksz domain is known to have broken this 
symmetry. However, this is similar to the RBC case 0] 
where the fluctuations continue to contribute a ring in 
Fourier space even above onset where convection rolls se- 
lect a direction in the plane. 

The radial structure of 5(k, Vo) seen in Fig. ^ is due 
to the optical transfer function T(k) of the shadow- 
graph, [lfj In Fig.Et we show the background-corrected 
azimuthal average SS(k, Vo) — S(k, Vq) — S(k, 0) of 
5(k, Vb). Figure^ gives T(k) calculated from Eq.|3]with 
the parameters of our optical system. |20| It contains the 
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FIG. 2: (a) and (b): Shadowgraph structure-factor 8S(k, Vo) 
as a function of the wave number k. (a): camera exposure 
time r = 0.1 s and (from top to bottom) Vo = 3.341,3.256, 
and 3.170 Volt, (b): V = 3.499 Volt and r = 0.500 (solid cir- 
cles), 4.000 (open circles), and 8.000 s (solid squares). Plusses 
near open circles: r = 4.000 s, / = 1000 Hz. (c): optical 
transfer function Eq. |!|] Dotted line in (a): fit of Eq. Qto 
the data for V = 3.341 Volt. Solid line: B 2 S n (k) (Eq. 
corresponding to that fit. 



oscillations of SS(k) observed in the data. 
We fitted the equations 
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to the data for SS(k,Vo), least-squares adjusting B and 
£. Thus, we used the Swift-Hohenbcrg form Eq. [5] for 
the structure factor S n of the refractive index and the 
optical transfer function [jj [Jjj in the weak-diffraction 
limit Eq. In Eq. z\ is the effective optical distance 
of the shadowgraph, = 2/2/ where I is the pinhole 
size of the light source and / is the focal length of the 
collimating lens, K is the wave number of the light, B = 
2dK, and Ji is the first-order Bessel function of the first 
kind.0,H3| As an example we show the fit to the data for 
V = 3.341 Volt in Fig. 01 as a dotted line. The solid line 
in Fig. gives the corresponding result for B 2 S n from 
Eq.0 

Results for B are shown in Fig. [3] as open squares. 
They vary smoothly through the transition, continue to 
rise beyond it, and finally have a maximum about 10% 
above Vf- A very different picture emerges from the data 
for £ that are shown in the insert in Fig. [3J We see that 
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FIG. 3: The amplitude B (open squares) and power P = 
7T 2 (solid circles) of the refractive-index structure factor 
S n , as a function of the voltage amplitude Vo. The insert gives 
the correlation length £ [inverse width at half height of S n (k)] 
as a function of Vo. Dotted lines: Vf = 3.22 Volt. 




FIG. 4: The relaxation rate T(k) as a function of the square 
of the wave number k. Solid circles: Vo = 3.50 Volt. Open 
circles: V = 3.25 Volt. Solid squares: V = 3.00 Volt. The 
lines are fits of T = To + T~2k 2 + Yik* to the data. Insert: 
example of SS(k, Vo, r) on a logarithmic scale as a function of 
r on a linear scale for Vo = 3.50 Volt, k — 1.327, and / = 100 
Hz (solid circles) and 1000 Hz (open circles). Solid line: fit of 
Eq.gjto the 100 Hz data. 



£ has a sharp cusp at Vf = 3.22 Volt, indicating that 
the transition is sharp roughly at the one percent level. 
However, there is no divergence of £ as would be expected 
for a second-order phase transiton. From £, B, and Eq.|2 
we can compute the total refractive-index fluctuation- 
power P = 2ir r o °° kS n {k)dk = ir 2 B/2£ 2 . This is shown 
as solid circles in Fig. |21 The Vo-dependence of P is 
similar to that of B\ no singularity at Vf is apparent 
from the data. Beyond Vo ~ 3.6 P decreases again. 

In order to learn about the dynamics of the fluctua- 
tions, we took images with several camera exposure times 
t at each of several values of Vb . Larger values of r lead to 
more averaging of the fluctuations and thus to a smaller 
5(fc, Vb,r). This is illustrated by the results shown in 
Fig. We are not aware of a prediction of the form of 
the time-correlation function C(k,St) of the fluctuations 
in a thin nematic sample and for Vo > 0, such as was done 
for RBC 0|- Thus we assume that C(k,St) is well ap- 
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FIG. 5: The relaxation rate To (solid circles: f = 100 Hz; 
solid squares: f = 1000 Hz) and the inverse square of the 
correlation length l/£ 2 (open circles, f=100 Hz) as a fuction 
of V 2 . Both have a cusp but remain non-zero at Vf The lines 
are straight-line fit to the data near the transition. 

proximated by a simple exponential decay exp[— T(k)Si\. 
One can show that this leads to 

5S {k,V ,T = SS (k,V o ,0) tvTITI ( 4 

r(fe, v )T 

for the average over a time interval r. At each of several 
Vo, we fitted Eq. 0] to measurements of SS(k, Vo, t) at 
each of many values of k to give SS(k, Vo, 0) and T(k, Vo). 
An example is shown in the insert of Fig. Although 
SS(k, Vo, 0) is influenced by T{k), we note that T(k, Vo) 
is not. In Fig. 0]we show some of the results for T as 
a function of k 2 . The data can be represented well by 
r = r + T 2 k 2 + r 4 fc 4 (solid lines in the figure). Results 
for To are given in Fig. [5] One sees that To has a sharp 
down-pointing cusp at Vf (vertical dotted line); but con- 
trary to predictions for a second-order phase transition 
Tq remains finite. Figure [5] also shows the results for £ 
in the form of 0.05/£ 2 as a function of V 2 . For a second- 
order transition and near Vf one would expect to be 
a linear function of V 2 and to vanish at Vf- However, 
£ -2 behaves very similarly to To and remains finite. The 
solid and dashed lines are straight-line fits to the data 
near Vp, separately above and below the transition. The 
fits for V 2 < Vp extrapolate approximately to To = 
and l/£ 2 =0 at V 2 ~ 13 (V c ~ 3.6), consistent with a 
point of absolute instability at V c that, for a first-order 
transition, can not be reached in the presence of strong 
noise. Interestingly, B and P (Fig.|3|) reach their maxima 
near V c . 

We presented data for both the statics and the dynam- 
ics of the bend Freedericksz transition of a homeotrop- 
ically aligned NLC that are inconsistent with a second- 
order phase transition and suggestive of a first-order tran- 
sition in the presence of strong noise. This disagrees with 



our theoretical understanding of this system that pre- 
dicts a second-order transition at the mean-field level. 
A fluctuation-induced first-order transition 0] should 
manifest itself only in a very small region, within 0.1 
percent or so of the mean-field transition, whereas the 
data suggest a critical region of width e c = O(0.1). 

This work was supported by the National Science 
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